METHOD AND APPARATUS FOR MEASURING ELASTIC 
CHARACTERISTICS OF A MEDIUM BY USING ULTRASOUND IMAGES 

FIELD OF THE INVENTION 

[0001] The present invention relates generally to the area of ultrasound signal 
processing, and more particularly, to a method and apparatus for measuring the elastic 
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characteristics of a medium by estimating variations of speckle patterns in ultrasounds 
images. 

BACKGROUND OF THE INVENTION 

[0002] Ultrasound imaging apparatuses are widely used, for medical diagnosis 
purposes, to obtain an image showing the internal organs of a human body. The 
ultrasound imaging apparatus transmits ultrasound signals toward a target object to be 
diagnosed, and utilizes the echo signals reflected from the target object to obtain an 
image showing the desired target object. One of the obstacles encountered with a 
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conventional ultrasound imaging apparatus in visualizing pathological tissues, such as 
tumors and cancerous tissues, is that if reflectivity between the pathological tissues and 
the neighboring tissues are not significantly different, then the conventional ultrasound 
imaging apparatus cannot provide a clear enough image of the pathological tissues for 
diagnosis. Conventionally, to overcome the above obstacle, a method of visualizing 
the elastic characteristics of an entire medium, including the pathological tissue and its 
neighboring area, is used to discriminate between the pathological tissue and the 
normal tissue. The method exploits the pathological tissue's property of having a 
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greater stiffness than the neighboring normal tissues. 
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[0003] Two kinds of methods for measuring the elastic characteristics of a medium are 
known in the art: one of the methods is to externally apply pressure to a medium of 
interest, measure the physical displacement of the medium, and calculate the local 



displacement from the pressure applied to the medium; and another exploits the 
difference in wave transmission properties between a stiff medium and a soft medium, 
and consists of the steps of applying low frequency vibrational forces to a medium and 
visualizing movement of the applied low frequency wave. However, all these 
methods have drawbacks as explained below. 

[0004] In the first method, a relative displacement must be estimated between the echo 
signals received in response to transmission of an ultrasound signal to the medium 
after applying pressure thereto and before applying pressure to the medium, in order to 
diagnose if the medium is pathological tissue. Preferably, the pressure applied to the 
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medium is increased since accurate diagnosis requires information on displacements in 
response to a wide range of pressures. However, if the pressure is increased, the 
correlation between the echo signals before and after applying pressure decreases, such 
that correctly identifying the displacement is difficult. As a result, discriminating 
between normal and pathological tissues is not viable. Further, in both methods, 
calculating waveform displacements or phase changes based on RF or envelope data, 
relating to the transmitted and received ultrasound signals is necessary. Since these 
methods need to compute signal correlations, they require complex hardware having 
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high-speed signal processing capability, and thus the costs associated with these 
methods is high. 

[0005] Accordingly, a need for a new and improved method and apparatus for 
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measuring the elastic characteristics of a medium exists in the art. 
SUMMARY OF THE INVENTION 

[0006] It is therefore an object of the present invention to provide a method and 
apparatus for measuring the elastic characteristics of a medium by estimating 
variations of speckle patterns in ultrasound images. 
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[0007] It is another object of the present invention to provide a method and apparatus 
of measuring the elastic characteristics of a medium, which can minimize the required 
number of frame memories and the required amount of calculation. 
[0008] According to one aspect of the present invention, a method for measuring 
elastic characteristics of a medium is provided, which comprises the steps of applying 
vibrations to the medium; acquiring a plurality of ultrasound image frames of the 
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medium; estimating a variation in brightness of a speckle pattern over the plurality of 
ultrasound image frames; and measuring the elastic characteristics of the medium 
based on the estimated brightness variation. 

[0009] According to another aspect of the present invention, an apparatus for 
measuring elastic characteristics of a medium is provided, which comprises a vibrator 
for applying vibrations to the medium; transducers for acquiring a plurality of 
ultrasound image frames of the medium; means for estimating a variation in brightness 
of a speckle pattern over the plurality of ultrasound image frames; and means for 
measuring the elastic characteristics of the medium based on the estimated brightness 
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vibration. 

BRIEF DESCRIPTION OF DRAWINGS 

[0010] The features of the present invention, which are believed to be novel, are set 
forth with particularity in the appended claims. The present invention, both as to its 
organization and manner of operation, together with further objects and advantages 
thereof, may best be understood with reference to the following description, taken in 
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conjunction with the accompanying drawings. 

Fig. 1 shows a distribution of reflectors within a sample volume of ultrasonic 
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beams transmitted from a transducer array. 

Fig. 2 is a view to illustrate a variation of a distribution of reflectors within a 



tissue in response to an applied vibration. 

Fig. 3 illustrates pixels situated in the same position in consecutive B-mode 
ultrasound images. ' " . 
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Fig. 4 is a schematic view showing a low frequency vibrator brought into 
contact with a human body, in accordance with the present invention. 
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Fig. 5 is a circuit diagram for calculating variances of brightness of pixels over 
a plurality of image frames, in accordance with the present invention. 
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DETAILED DESCRIPTION OF THE PRESENT INVENTION 

„..- ■ • • . 

[0011] Ultrasound images generally include a particular noise pattern called a "speckle 
pattern." For example, a B-mode ultrasound image of a soft tissue within a human 
body has a plurality of granular speckles. The speckles are generated by the coherent 
nature of the ultrasound imaging system and are a detriment to obtaining high- 
resolution ultrasound images. More specifically, the speckles are generated because 
the signals reflected from the reflectors distributed within ultrasound beams are non- 
uniformly superimposed. Therefore, as shown in Fig. 1, if the distribution of 
reflectors within a sample volume of ultrasound beams from the array transducers 
varies, then the speckle pattern varies accordingly. Therefore, if a vibrational force is 
applied to soft tissue of the human body to vibrate the tissue, the speckle pattern in the 
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ultrasound image changes. Thus, the ultrasound image brightness changes. 
[0012] As mentioned before, the stiffness of pathological tissue is greater than that of 
normal soft tissue. Further, in response to an applied vibration, the change in 
distribution of reflectors within the tissue is inversely proportional to the stiffness of 
the tissue, as shown in Fig. 2. In other words, for stiff tissue, an applied vibration 
causes the entire tissue, including the tissue of interest, to move as a whole, such that 
the speckle pattern does not change much. On the other hand, for soft tissue, the 
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applied vibration significantly changes the positions of the reflectors within the tissue. 
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As a result, the change in the speckle pattern, in response to an applied vibration is 
proportional to the softness of the tissue. Based on this principle, the stiffness of the 
tissue may be measured. 

[0013] The change in brightness of a speckle can be estimated by observing the 
variation of pixels situated at the same positions along the time axis in a series of B- 
mode ultrasound images. As an indicator of a change in brightness, a variance of 
brightness may be used. Alternatively, a standard deviation, which is the square root 
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of a variance, may also be used. According to the present invention, variances of 
brightness for all the pixels in an image frame are calculated, and their values are 
displayed in real time for the measurement of the elastic characteristics of a medium. 
The user may arbitrarily determine the number of image frames for calculating 
variances. In order to obtain a more accurate estimate of the elastic characteristics of 
a medium, the magnitude and frequency of an applied vibrational force may be 
adjusted to a find a resonant frequency and magnitude. 

[0014] The vibrator frequency and amplitude may be adjusted is needed. Preferably, 
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the frequency of the vibrator is within several hundred Hertz. Fig. 4 shows a 
schematic view where a low frequency vibrator is brought into contact with a human 
body. Referring to Fig. 4, an operator diagnosing a patient applies vibrator 15 to the 
patient's body. As shown, the vibrations are applied in a non-invasive manner by 
bringing vibrator 15 into contact with the skin of the human body and keeping 
transducer 10 down toward the human body. The tissue vibrates according to the 
force applied by vibrator 15, and as a result, scatterers consisting of the tissue move so 
that a plurality of different images frames as shown in Fig. 3 are obtained. Thus, the 
brightness of a certain pixel may continuously changes throughout the image frames. 
[0015] Alternatively, an invasive method may be used to transfer the vibrational force 
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to a medium to more precisely observe specific tissue of the human body. A needle is 
inserted into the medium and vibrations are transferred to the medium via the tip of the 
needle. Although this invasive method may be somewhat burdensome to the patient 
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as the needle has to be inserted into the body of the patient ■ the advantage is that 
vibrations may be precisely transferred to a desired location in the body. 
[0016] The speckle brightness change due to an applied vibration may be determined 
by calculating the brightness change of each pixel. As a measure of the brightness 
change, the brightness standard deviation at each pixel of an entire B mode image 
frame may be used. A mean m n (/, j) and a standard deviation a n (/,.-/) of a certain 
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number, i.e., M, of consecutive image frames, starting from the (/z-M+1) frame and 
including up to the n th frame, at individual pixel points (i, j) may be expressed by 
Equations (1) and (2), respectively: 
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Eq- (i) 

M k-n-M + \ 



1 n 

a ^^ = "\ll7 %[ x ^J)- m n(jJ)\ Eq.(2) 



wherein x k (i,J) is the brightness at pixel (i, j) of the k th frame, which corresponds to 
each pixel value displayed on a monitor; and i and j are horizontal and vertical pixel 
indices, respectively. 

[0017] In calculating brightness variances for pixels in the image frame?, N is first 
assumed to be the current frame number, N-l is the previous frame number, X N (i, j) 
represents a brightness value of the (i, j) th pixel in the n th frame, and m N (i, j) and m N .!(i, 
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j) represent an average brightness for the (i, j) th pixel over the N number of frames as 
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counted from the first frame and an average brightness for the (i, j) pixel over the N-l 

number of frames as counted from the first frame, respectively. cr N (i, j) and 

&n-\ ('>/) are a l so assumed to represent brightness variances for the (i, j) pixel over 

th 

the N number of frames as counted from the first frame and for the (i, j) pixel over 
the N-l number of frames as counted from the first frame, respectively. Then, the 
average brightness for the (i, j) pixel, rrtN(i, j)> and the brightness variance for the (i, 
j) th pixel, a\ , may be expressed by Equations (3) and (4), respectively: 



m N (ij') = m NA (i,j)+^[x N (i,j)-m N ^(i,j)] 

A Eq. (3) 

t" N (hj)= — — — 

N 



°l (i> J) = (h j) + m 2 N _ x (i, j) + -J-[4 0', J) - Nml (i, j) - (ij) - m\_ x (i, j)} 

N 
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Eq.(4). 



Equations (3) and (4) are recursive forms of Equations (1) and (2), respectively. 
[0018] According to the above recursive equations!, the previous average and variance 
values are used to calculate the current average and variance values. The initial 
values for the average and variance are set based on the first inputted data pursuant to 
the general principle used in recursive types of equations. Therefore, the average and 
variance in the first frame may be represented by Equation (5): 



Eq. (5). 



[0019] That is, the initial value for the average, mi(i, j), and initial value for the 
variance, a] (ij) , are equal to Xi(i, j) and 0, respectively. The next values for the 
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average and variance are recursively calculated based on the next input X2(i, j). If the 
average and variance of brightness are recursively calculated according to above 
Equations (3) and (4), images showing variance values of brightness for all the pixels 
may be displayed at the same rate as the frame rate of the input images since the 
variance values are calculated upon acquisition of a new frame. 

[0020] Referring to Fig. 5, there is shown a circuit diagram for calculating variances of 
brightness for pixels over a plurality of image frames, according to the above equations. 
As shown in Fig. 5, a circuit 600 includes a plurality of adders 601, a plurality of 1/N 
dividers 602, a plurality of frame delay elements 603, a plurality of squarers 604, and 
N-fold multiplier 605. For frame delay element 603, a frame memory may be used. 
The first frame delay element is used for storing the average taken up to the (N-l) th 
frame, mN-i(i, j), and the second frame delay element is used for storing the variance 
taken up to the (N-l) th frame, a 2 N _ } (ij) . According to the structure of the circuit for 
calculating the variance as shown in Fig. 5, only the average and variance for each 
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pixel in the image frame that is immediately prior to the current image frame are stored. 
The circuit of Fig. 5 may be implemented by either hardware or software, and is 
preferably mounted at the rear end of a conventional ultrasound imaging system along 
with the vibrator. 

[0021] In order to further reduce the amount of calculations, an absolute difference 
between maximum and minimum values of brightness of a pixel over a predetermined 
frame period may be used instead of the variance, which has to be calculated at every 
frame timing. Assuming, for example, 1 5 image frames Xi(i, j)-Xi 5 (i, j), maximum 
and minimum values of brightness of every pixel over the 10 image frames, i.e., maxj. 
io(i, j) and mini_io(i, j) ? may be calculated. Thereafter, a difference between the 
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maximum and minimum values, i.e., diff x _ l0 (i,j) = max 1 _ 10 (z 5 y)-min 1 _ 10 (/" 5 y)| may 
be calculated. Performing this type of calculation for every pixel results in a 
difference image diff\.\o, indicating the difference at every pixel position. Likewise, 
if the calculation is performed upon the image frames X2(i, j)-Xn(i, j), a difference 
image diff2-\ i may be obtained. Repeating the calculation finally results in a 
difference image differs- Repeating the same calculation for the next image frames 
Xi6(i> j), Xiy(i, j), etc., provides a difference image updated at every frame timing. 
According to this embodiment, an image similar to the variance image may be 
obtained with a reduced amount of calculation as compared to the method illustrated in 
Fig. 5. / 

[0022] As explained above, the reason why the variance of brightness over the image 
frames may be used as an indicator of the variation of the speckle pattern is that the 
variance of brightness for a certain pixel is proportional to the displacement of the 
scatterer located at the pixel position, and the displacement is proportional to the 
stiffness of the scatterer. Since pathological tissue has an increased stiffness as 
compared to normal soft tissue, the speckle pattern for pathological tissue does not 
vary much because the variances of brightness as calculated by using the above 
equations and the circuit shown in Fig. 5 do not vary significantly over several image 
frames. On the other hand, for normal soft tissue, the scatterers within the tissue 
move to a great extent under an applied vibration, and thus the variances of brightness 
calculated by using the circuit shown in Fig. 5 vary significantly over several image 
frames. Based on this principle, the stiffness of the tissue may be measured. For a 
more accurate diagnosis, the standard value for the stiffness of certain tissue must be 
provided so that the degree of deviation from the standard value may be used as a 
measure for the diagnosis. 

[0023] As described in the foregoing, the elastic characteristics of a medium may be 



measured by calculating variances of brightness for every pixel over consecutive B- 
mode ultrasound images. Therefore, the operator who is diagnosing the patient may 
readily discern any pathological tissue within the patient's body. According to the 
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method and apparatus of the present invention, only acquired ultrasound images are 
used for the measurement, and thus software programs may be employed for such 
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purpose. Further, since the method and apparatus of the present invention employs a 
recursive type of calculation, the size of the required hardware such as frame memory, 
amount of calculation, and the CPU load may be greatly reduced. 
[0024] While particular embodiments of the preseint invention have been shown and 
described, those skilled in the art will appreciate that changes and modifications may 
be made without departing from this invention in its broader aspects and, therefore, the 
aim in the appended claims is to cover all such changes and modifications, as fall 
within the true spirit and scope of this invention. 
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